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Abstract 

Four-fermion processes with a particle lost in the beam pipe are studied at LEP to 
perform precision tests of the electroweak theory. Leading higher-order qed correc- 
tions to such processes are analyzed within the framework of the Structure Functions 
(sf) approach. The energy scale entering the qed sf is determined by inspection of 
the soft and collinear limit of the 0(a) radiative corrections to the four-fermion final 
states, paying particular attention to the process of single- W production. Numerical 
predictions are shown in realistic situations for lep experiments and compared with 
existing results. A Monte Carlo event generator, including exact tree-level matrix 
elements, vacuum polarization, higher-order leading qed corrections and anomalous 
trilinear gauge couplings, is presented. 
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1 Introduction 



Four-fermion final states are of special interest for the physics programme of 
LEP2 and future high-energy electron-positron colliders, being entangled with 
electroweak gauge boson production and decay [1]. In particular, the process 
considered in the present paper, i.e. 

e + e" -> e~(e+)z/ e (z/ e ) q'q (1) 

is peculiar among all the possible four-fermion final states because the bulk 
of its cross section is due to two sub-processes, i.e. W^-boson pair production 
and decay 

e + e" -> W*W* -> 4 fermions (2) 

and the radiation of an almost on shell t-channel photon from the electron 
(positron), with subsequent production of a W^-boson and a neutrino 

e + e~ — > 7*e + (e~) — > W*v e {e~) — > 4 fermions (3) 

Despite, strictly speaking, both sub-processes (2,3) always occur simultane- 
ously and are indistinguishable, channel (2) dominates if the electron is emit- 
ted at large angle, whereas channel (3) dominates if the electron is emitted 
in the very forward region, because of the presence of a quasi-real t-channel 
photon. 

In this paper the process (3) will be addressed, by restricting the analysis 
to the kinematical range of forwardly emitted electrons. This signal is usually 
referred to as single- production, since only the two final jets are detected [2]. 

The importance of this process has been emphasized since long time. In the 
lep2 energy range it is fundamental in order to study the self-interaction 
of the gauge bosons, together with the process (2), whereas in the energy 
regime of future colliders at the TeV scale it becomes the dominant elec- 
troweak process. In refs. [3,4] cross sections and distributions were calculated 
in the approximation of real W^-boson production, either by [3] studying the 
reaction e + e~ — > e^u e W + , or by [4,5] employing the Weizsacker- Williams [6] 
equivalent-photon approximation for the t-channel photon. In refs. [3-5] it was 
pointed out the relevance of this process for the study of trilinear gauge boson 
couplings and some assessment of the sensitivity has been given. The first full 
four-fermion calculation, including the crucial effect of fermion masses, has 
been presented in ref. [7], where the lep2 sensitivity to anomalous gauge cou- 
plings has been studied. Since then, other complete four-fermion calculations 
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of the single-VT process have appeared in the literature and implemented in 
computational tools for data analysis [7-11]. In most of these calculations the 
effect of fermion masses is exactly accounted for in the dynamics and kinemat- 
ics for the whole four-fermion phase-space [7,8,10,11], while in the approach 
of ref. [9] the Weizsacker- Williams approximation is employed in the very for- 
ward, collinear region and matched with a massless four-fermion computation 
outside it. An up-to-date inventory of the present theoretical status is under 
preparation by the four-fermion working group of the LEP2 MC workshop at 
CERN [12]. 

Measurements of the single-VT cross section and the corresponding bounds on 
anomalous gauge couplings have been recently reported by the lep collabo- 
rations [2]. Because the foreseen accuracy of final lep2 data is of the order of 
1-2% [12], accurate theoretical predictions for cross section and distributions 
are required. 

The calculation of the cross section for single-Vr processes poses several non- 
trivial theoretical problems [12]. For a realistic account of gauge bosons prop- 
erties it is mandatory to include the gauge boson width in the propagator. In 
general this mixes a fixed order calculation with an all order resummation of a 
class of Feynman diagrams and introduces a violation of the Ward identities of 
the theory. This issue is of special importance here since, due to the t-channel 
photon exchange, even a tiny violation of QED Ward identities is enhanced 
by a factor of s/m 2 e . This is indeed the case if a running width is used in the 
calculation. This problem has been extensively studied [13-15], and several 
options to address it have been explored. The most theoretically appealing 
procedure is the fermion loop scheme [14,15], which preserves both U(l) and 
577(2) Ward identities. Recently, this scheme has been generalized to the case 
of massive external fermions, both in its minimal version, which considers the 
imaginary parts of the fermionic loops (ifl) [8], and in its full realization with 
real and imaginary parts [16,17]. In particular, in ref. [8] a detailed numerical 
investigation has been performed, showing no significant difference between 
the ifl and the fixed width scheme, even in the region most sensible to U(l) 
gauge invariance. For this reason, the fixed width scheme is adopted in the 
present calculation. 

Another delicate issue is the so called resolved-photon component of the cross 
section. The quasi-real t-channel photon can split into a pair of almost mass- 
less quarks, leading to a situation where the partonic picture of hadrons is 
inadequate and both perturbative and non-perturbative QCD corrections be- 
come relevant. This issue is widely discussed in the literature [18], where the 
standard approach to this problem is also described and to which the reader 
is referred for details. However, for single-W^ like events the resolved-photon 
component does not constitute a severe limitation: once a hard qq invariant 
mass cut is imposed, as done in realistic data analysis, the bulk of the signal is 
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kept, whereas the resolved photon contribution becomes almost negligible [9]. 

A further relevant issue is given by radiative corrections due to photon radia- 
tion. Because exact 0(a) electroweak corrections to single- W production are 
still unknown, in most of the theoretical and experimental studies presented in- 
sofar only the large contribution of initial-state radiation (isr) has been taken 
into account, generally by using collinear structure function (hereafter denoted 
as Sf) and assuming s = 4E 2 as the proper scale for qed radiation. Due to the 
dominance of the quasi-real t-channel photon exchange, this can be expected 
not to be a suitable choice in the present case. On the other hand, it has been 
recently proposed to correct only the s-channel contributions to the single-14 7 
signature, fixing the radiation scale in the usual manner, and to neglect the 
photonic corrections to the t-channel contributions [10]. Following previous 
investigations [19-23] of the pattern of photonic radiation in qed and elec- 
troweak processes, some theoretical arguments to determine the appropriate 
energy scale entering the SF are presented and compared with existing results. 
The analysis here described elucidates the theoretical details and provides fur- 
ther numerical results of a contribution by the authors [24] to the activity of 
the four-fermion working group of the LEP2 MC workshop at cern [12]. Ideas 
similar to those adopted in the present work have recently appeared in ref. [25] 
and there applied to the two-photon process e + e~ — > e + e~fi + fi~. 

The paper is organized as follows. After a short review of the SF approach 
to leading log (LL) QED corrections in Sect. 2, Sect. 3 collects the analytical 
results valid for soft and collinear corrections to a generic scattering process. 
By comparing the results of Sect. 2 and Sect. 3, the radiation scales for single- 
W production are determined in Sect. 4. Sect. 5 deals with the problem of 
taking into account the effect of the photon vacuum polarization in the single- 
W process, while Sect. 6 shows the numerical results of the present study 
obtained with a Monte Carlo (hereafter mc) code for the single-VF signature, 
including also the effect of anomalous trilinear gauge couplings. Conclusions 
and prospects are given in Sect. 7. 



2 Structure Function approach to photon radiation 

Since in high-energy processes the corrections due to the emission of soft and 
collinear radiation are quite large, the LL contribution must be calculated 
at every perturbative order. A common technique to achieve this goal is the 
QED structure function approach [26], which consists in convoluting the hard- 
scattering cross section with appropriate "parton" densities. As well known, 
these convolution factors, i.e. the QED structure functions, include, by con- 
struction, both the real and virtual part of the photon correction, in order 
to ensure the cancellation of the infrared singularities. If a generic Born-level 
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prediction dao is considered, the cross section da including LL qed radiative 
corrections is obtained, by virtue of factorization theorems, according to the 
following general formula [26] 



da — Y[ / dXiD(A 2 ,Xi) da { 



o 



(4) 



where 1 Obi SX6 the energy fractions carried away by the radiated photons 
from the i-th leg, A is the characteristic scale of the SF D(A 2 , Xj), whose evo- 
lution is driven by the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) 
equation [27] and is dependent on A. It is worth noticing that the choice of 
the scale A is not dictated by general arguments and it is therefore rather 
arbitrary. 

Equation (4) can be rewritten by stressing the possibility of different scales 
for each SF as follows 



In particular, if the integrated hard-scattering cross section is a smooth func- 
tion of the centre of mass (cm.) energy, once the integrations over the energy 
fractions Xi are performed in the soft-photon approximation, the 0(a) double- 
log expansion of eq. (5) can be written as follows 



where AE/E is the maximum energy for undetected photons, to be identified 
with finite energy resolution of the photon detector, and L(Af) = log(A^/m 2 ) 
is the collinear logarithm. 

Since the functional form of the QED SF is accurately known [26], the main 
problem in evaluating eq. (5) is to fix the process scales Aj. A generally adopted 
attitude is given by the comparison of eq. (6) with a perturbative calculation, 
which can be performed within any approximation, provided it reproduces the 
correct double-log contribution of the 0(a) correction. This issue is addressed 
in the next Section. 



3 Analytical results 

The double-log contribution to photon radiation traces back to soft and collinear 
bremsstrahlung and its virtual counterpart [28], and, in the case of a calori- 
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metric measurement of the energy of the final-state (fs) particles, to hard 
radiation collinear to the FS particles themselves [29-32]. At each perturba- 
tive order, the leading contribution can be expanded in terms of infrared and 
collinear logs. For example, when the photons are emitted from the initial-state 
(is) particles only, such an expansion can be arranged in terms of double-log 
contributions of the form a n l n L^, where I = \og(AE/E) is the infrared log 
and L s = log(s/m 2 ) is the collinear log. This is the reason why A 2 = s is the 
"natural" energy scale to be used for SF in the presence of ISR only. When 
also FS radiation is considered, the collinear log, L, is in general modified by 
additional factors coming from the angular integration over the photon vari- 
ables. A typical example is given by the radiation emitted from one leg in 
the t-channel QED contribution to Bhabha scattering [19]. In the soft-photon 
approximation the radiation cones, one from the IS electron and one from the 
FS electron, have a half-opening which is determined by the angle between the 
emitting particles, because of a destructive interference. As a consequence, 
the energy scale s, which appears in L s , transforms into \t\ ~ s(l — cos8), 
where is the electron scattering angle, and, therefore, L s — > L t . Hence the 
perturbative expansion contains collinear logs which are modified because of 
the angular ordering introduced by the radiation cones. In the presence of 
large scattering angles, for which \t\ ~ s, the above modification is numeri- 
cally small, but it becomes more and more important in the forward angular 
range, which is the dynamically favourite region by t-channel Bhabha scatter- 
ing and where t C s. The net result is a numerically significant depletion of 
QED radiation effects just in the most important part of the hard-scattering t- 
channel dynamics. Actually, when using SF to evaluate QED LL corrections to 
small-angle Bhabha scattering, the energy scale A 2 = \t\ is employed in all phe- 
nomenological applications [19]. More in general, in order to take into account 
dominant initial-final-state interference effects in addition to initial- and final- 
state leading terms, s and t QED contributions to Bhabha scattering can be 
corrected in terms of a unique combination of Mandelstam invariants given by 
st/u, as discussed in refs. [20,21]. Therefore, the energy scale A 2 = st/u turns 
out to be a suitable choice for the evaluation in terms of SF of LL corrections 
to QED Bhabha scattering, as demonstrated, in comparison with the exact 
0(a) calculation, in ref. [21]. Similar arguments for an appropriate choice of 
the energy scale for QED radiation, based on the inspection of the soft and 
collinear limit of the 0(a) correction, have been also advocated in ref. [22] for 
the reaction e + e~ — > W + W~ and, very recently, in ref. [25] for the process 
e + e~ — > e + e~/j + fi~. Comparisons performed in refs. [22,25] with available ex- 
act 0(a) calculations explicitly exhibit the validity of such a strategy, which is 
therefore pursued in the present analysis. As already remarked, the result for 
LL corrections in the presence of a calorimetric detection of FS particles must 
include the contribution of photons which, regardless of their energy, can not 
be discriminated from closely collinear fermions, as a consequence of the finite 
angular resolution of the calorimeters. The role of such hard photons collinear 
to the FS particles becomes therefore unavoidable in the case of a calorimetric 
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measurement of the energy of the FS particles, as discussed in the following. 



3.1 Soft-photon contribution 



In this Section the contribution of photons, too soft to be detected in the 
calorimeter, will be computed for a generic process with n ingoing-legs. Q The 
following approximations are understood 

« » * (7) 

3> m? m 2 A 



where is the momentum of a particle of mass emitting a real photon of 
momentum k, and Sy = (<& + g,,) 2 is the invariant mass of the pair ij. 

By following the standard derivation of the eikonal factors due to soft brems- 
strahlung and by generalizing it to particles with different masses and charges, 
the differential cross section, dressed by soft-photon emission, can be cast into 
the following factorized form [28] 

du) 2a J\ _ _ . s,^ . . 

acTsoft = da > QiQj log (8) 



where a; is the photon energy and Qi is the charge of the z-th particle. 

It is worth noticing that in the limit <C ?^ 2 , m 2 , provided the first in- 
equality in eq. (7) still holds, the logarithmic behaviour present in eq. (8) 
disappears, leaving a power law which can be simply obtained by means of 
the substitution [28] 

log^^i^ (9) 

miTJij 6 miTJij 



Notice that, since the goal is to determine the scale entering the SF, only 
the contribution of real photons is explicitly calculated, because the virtual 
corrections, in order to preserve the cancellations of infrared singularities, must 
share the same leading collinear structure of the real part itself. 

By including the virtual part needed to cancel the infrared singularity and 
integrating eq. (8) over the photon energy u in the soft region < uj < AE, 

1 This choice fixes our conventions. Outgoing particles will appear as ingoing ones 
with momentum and charge according to crossing symmetry. 
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one gets 



AE2a " Sij 
d<r s +v = ao-o log -j^-— Z^QiQj lo S ; 



£ vr 



mirrij 



(10) 



3.2 i/ard radiation collinear to the final-state particles 



In the case of a calorimetric set-up, which is the realistic situation for single- W 
production at LEP, photons collinear to the detected FS particles can not be 
distinguished from the the emitting particles themselves, independently of the 
photon energy. Therefore, in order to obtain the correct structure of double- 
log corrections for such an event selection, the effect due to the emission of 
unresolved hard radiation collinear to the FS particles must be taken into 
account in addition to soft+virtual corrections. 



To this end, let us re-consider the previous process with n ingoing legs and 
think m of them to be changed into outgoing legs at the end of the calculation 
(see the previous footnote). Then, the contribution of photons collinear to the 
FS particles can be cast into a gauge invariant form as follows [30-32] 

du2a ™ 2 EiS 
o(7 hard = da 2^ Qi lo S ( 1 1 ) 

UJ 7T , TUi 



where Ej, is the energy of the i-th particle, and 5 is the half-opening angle of 
the calorimetric resolution. 



By integrating eq. (11) over the photon energy u in the range AE < uo < E, 
the integrated correction due to hard photons collinear to the FS particles is 
given by 

E 2a A „ 2 , EiS 

dffhard = da log — — — }^ Q { log (12) 

I\E 7r ^ mi 



3.3 The master formula 



Equations (10) and (12) give the leading double-log contribution which must 
be compared to (6), the 0(a) perturbative expansion of eq. (5), in order to 
fix the process scales A,. Summing the contributions of eq. (10) and eq. (12), 
the analytical cross section is in conclusion given by 
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lot A PI [ n E- 
da s+Y + da h&rd = da — log — - < V Q\ log — + (13) 

n m 

- £g i g,-io g 2(i- CiJ -)-i;^ lo g <J 

i>j i 



where is the cosine of the angle between particles i and j. 

Three different kinds of logarithms occur in eq. (13). The first term contains 
the mass and energy logarithms of the IS particles only, since, as expected, the 
energies and the masses of the FS particles disappeared, in agreement with the 
kln theorem [33]. The second term includes angular terms due to radiation 
interference, while the third one comes from the requirement of calorimetric 
measurement. 

These terms must be compared with the collinear logarithms of eq. (6) in order 
to fix the scales Aj of the SF. In the following Section this task is accomplished 
in detail for the single-IF process. 



4 Fixing the radiation scales in the single-VT process 



Let us consider, for definiteness, the process e + e~ — > e~uud with the FS elec- 
tron lost in the beam pipe (single- W process). In this event selection (here- 
after es) the leading contribution comes from 7*e + scattering with the virtual 
photon emitted from the electron line. The leading dynamics is given by the 
t-channel Feynman diagrams shown in Fig. 1. 



Fig. 1. The fusion and bremsstrahlung diagrams are the leading Feynman graphs 
for the single-VF signature. 

If a calorimetric measurement of the energies of the FS particles is performed, 
only the IS legs need to be corrected by the SF. Furthermore, since the electron 
is scattered in the very forward region, the interference between the electron 
line and the rest of the process is very small. This allows a natural sharing of 
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the logarithms coming from eq. (13) between the two IS SF associated to the 
colliding electron and positron. 



Hence the formula (13), when compared with eq. (6), translates into the two 
following scales (A_ refers to SF attached to the electron line, while A + to the 
SF attached to the positron line), 



where E is the beam energy, c_ is the cosine of the electron scattering angle, 
c u and cj are the cosine of the quark scattering angles with respect to the 
initial positron, and c u j is the cosine of the relative angle between the quarks. 

It is worth noticing that in the numerical implementation, whenever one of 
the two scales is less than a small cut-off (A^ ut _ off = Am 2 e) where m e is the 
electron mass), the radiation from the corresponding leg is switched off, in 
accordance with eq. (9). It was carefully tested that variations of the cut-off 
do not alter the numerical results. 

Owing to the presence of a resonant W boson, some modifications to the 
previous results may come from finite-width effects and from radiation deco- 
herence [34]. Finite- width corrections of the form of E 1 /Y\y arise when the 
unstable particle propagator goes off its mass-shell, but this is not the present 
case, since the multi-fermion final state can accommodate a resonant W. Ra- 
diation decoherence is present whenever a resonance occurs and its effect is 
to cancel the angular dependence from the scale. As a consequence the scale 
A + should be modified by dropping the angular interference factors in eq. (14) 
when the emitted photons have E 1 ~ Tyy. Yet in the present case the effect is 
tiny, since the effects due to angular interference for the scale A + are already 
small by themselves. 

It is also possible to make a naive ansatz for the radiation scales without 
a detailed calculation, by thinking of the graphs of Fig. 1 in terms of the 
Weizsacker- Williams approximation [6], i.e. in terms of a convolution of the 
process e + 7 — > u e W* with an equivalent photon spectrum plus a real electron 
line. This leads to assigning two different scales to the single-VF process: one 
scale for the electron current and one for the positron current. The former scale 
is the proper one for a t-channel process, e.g. t-channel Bhabha scattering, so 
it is simply |^*|, where \q^*\ is the squared momentum transfer in the eej* 
vertex. The latter is the sum of an s-channel electron exchange and a t-channel 
W exchange (see Fig. 1). Assuming that the t-channel dominates, its natural 
cut-off is given by the VF-boson mass, M w . Hence, the following ansatz follows 




(14) 



A 2 _ 



,naive ~~ l?7*l > ^+,naivc ~~ 



(15) 
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where Mw is the mass of the W boson. The comparison between the scales 
given by eq. (14) and these naive scales, which will be performed numerically 
in the following Section, provides a useful cross-check of the analytical results 
derived by inspection with the soft /collinear limit of the 0(a) correction. 

A discussion of other possible approaches to the treatment of photonic correc- 
tions to single-VT production can be found in the four-fermion working group 
report of the LEP2 MC Workshop [12]. 



5 The running of the electromagnetic coupling constant 

Besides the higher-order qed corrections discussed in the previous Sections, 
other large logarithmic contributions to the single-FY cross section arise from 
the running of the electromagnetic coupling constant a. Since in the case under 
study the dominant configurations come from the Feynman diagrams with 
an almost on-shell photon exchange, the appropriate scale for the evaluation 
of the electromagnetic coupling relative to the t-channel photon in the ee'j* 
vertex is the squared momentum transfer q 2 * defined above. 

However, because Gp, My/ and Mz are the typically adopted input parameters 
for electroweak processes at lep2, the electromagnetic coupling is fixed at 
tree-level to a high energy value as, for example, 

A /^^F^W S W -4.1, 2 -I na\ 

a Gp = 4V2 , with s = 1 - — — . (16) 

An M| 

On the other hand, the single-FY production is a q 2 * ~ dominated process 
and therefore the above high-energy evaluation of a, ac F , needs to be rescaled 
to its correct value at small momentum transfer. To this end, a gauge-invariant 
"reweighting" procedure can be adopted, by rescaling the differential cross 
section da/dt (t = q 2 *) in the following way 

da _^ a 2 jt) da 

dt ~^ a 2 Gp dt ' 1 ' 

where a(t) is the running coupling constant computed at virtuality q 2 *. 

A detailed analysis of the effect of the running couplings in single-H^ produc- 
tion has been recently performed within the massive fermion-loop scheme in 
ref. [17], where the couplings are automatically running in the calculation. As 
shown in ref. [17], the relative difference between the above reweighting pre- 
scription and the complete results of the fermion-loop scheme is at the l%-2% 
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level 0, and it is therefore in the expected range of theoretical uncertainty due 
to missing full one-loop electroweak corrections. 



6 Numerical results 

In this Section the MC code, developed to simulate the single-VT process, is 
described and a sample of numerical results obtained by means of it is shown 
and commented, with particular emphasis on the effects of higher-order qed 
corrections to single-W 7 production at LEP2 energies. 

6.1 The Monte Carlo code 

A MC program, named SWAP, was developed to calculate cross sections and 
differential distributions for the single-W 7 signature. 

As already emphasized, the main feature of this process is the fact that the t- 
channel photon of Fig. 1 becomes quasi-real. In the limit of massless fermions, 
the photon propagator becomes singular in the forward direction and the cross 
section develops a logarithmic singularity. Indeed, whenever the final electron 
is lost in the beam pipe, its mass becomes a natural cutoff for the very-forward 
singularities, compelling to build a massive matrix element and phase-space. 
The phase-space integration is performed in SWAP with the aid of a multi- 
channel importance sampling with stratification. The main peaking structures 
for the single-H^ process are given by the dynamics depicted by the fusion and 
bremsstrahlung graphs of Fig. 1. They are the resonant P^-boson invariant 
mass, treated with a Breit-Wigner weight, and the t-channel "singularity" of 
the quasi-real photon, treated with a l/|t| weight. Moreover, the program can 
deal with the singularities of the sub-leading t-channel CC20 diagrams shown 
in Fig. 2, by means of the multi-channel approach. 

The exact hard-scattering matrix element is computed by means of the ALPHA 
code [35] for the automatic evaluation of Born scattering amplitudes. Fermion 
masses are exactly taken into account and the fixed-width scheme is adopted 
as gauge-restoring approach, by taking the massive gauge boson propagator 
as follows: 



Actually, for the single-VK final state under examination here and for realistic 
event selections, the differences between the two procedures are confined below the 
1% level. 




(18) 



k 2 — M 2 + iTM 



12 



Fig. 2. The multiperipheral diagram is the main sub-leading Feynman graph for 
the single-PF signature. 

It is known [8,14,15] that this scheme preserves U(l) gauge invariance but 
still violates SU{2) Ward identities. However, at least in the unitary gauge 
employed here, it is indistinguishable from other fully gauge-invariant schemes 
[8,14,15]. 

The contribution of anomalous gauge couplings is also accounted for in SWAP. 
The anomalous gauge boson couplings Afc 7 , A 7 , 5z, Akz and A 7 are imple- 
mented in the ALPHA code according to the parameterization of refs. [38,39]. 
Photon radiation is implemented via SF formalism, according to the discussion 
of Sect. 4. It is worth noticing that, since the incoming electron/positron are 
required to be on-shell massive fermions, a naive four-momentum rescaling, 
due to photon emission, such as p± — xp± leads to potentially dangerous gauge 
violations, according to what previously discussed. Therefore, the rescaled in- 
coming four-momenta are implemented as p± = (x E, 0, 0, ±^Jx 2 E 2 — m 2 ), 
by interpreting x as the energy fraction after photon radiation, as motivated 
in ref. [21]. If required, p±/pl effects can be provided in the treatment of 
ISR, by means of either p_i_-dependent SF [36] or a QED Parton Shower al- 
gorithm [21,37]. The effect of vacuum polarization is taken into account as 
described by eq. (17), by including the contribution of leptons, heavy quarks 
and light quarks, the latter according to the parameterization of ref. [40] . The 
program supports realistic ES and it can be employed either as a cross-section 
calculator or as a event generator, with both weighted and unweighted gener- 
ation available. 

The technical precision of the event generator SWAP has already been carefully 
proved in ref. [12], by means of detailed tuned comparisons between the pre- 
dictions of independent codes. Perfect agreement was found, both at the level 
of integrated cross sections and distributions, also for purely leptonic final 
states. 
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Table 1 

The ES adopted for the calculations shown in the present paper for the signature 
e + e~ — > e~v e ud, according to ref. [12]. 



electron angular acceptance 
quarks angular acceptance 
calorimetric half-opening angle 


| cos 6» e | > 0.997 
1. no cut 
5.00° 


| cos 6» e | > 0.997 
2. | cos 6 q ,q\ < 0.95 
5.00° 


quark-antiquark invariant mass 


45.0 GeV 


45.0 GeV 



Table 2 

The adopted electroweak input parameters, according to ref. [12]. All other param- 
eters are calculated by means of the tree-level relations. 

G F = 1.16637 x 10- 5 GeV" 2 
M z = 91.1867 GeV 
M w = 80.35 GeV 
T z = 2.49471 GeV 



6.2 Discussion of the numerical results 

The numerical simulations are elaborated according to the ES reviewed in 
Tab. 1, with the electroweak input parameters shown in Tab. 2. 

In Figs. 3-4 the numerical impact of different choices of the A 2 -scale on the 
cross section of the single- W process e + e" — > e~vud in the lep2 energy range 
is shown. Since the energy scale A + of eq. (14) depends on the quark scattering 
angles, two different quark angular acceptances are considered, namely no cut 
(Fig. 3) and | cos$ U; j| < 0.95 (Fig. 4). The marker • represents the Born cross 
section, o represents the correction given by A 2 = s scale for both IS SF, 
represents the correction given by K\ = |g 2 »| scale for both IS SF, a the 
correction given by the scales of eq. (14), □ the correction given by the naive 
scales of eq. (15). It can be seen that neither the s scale, as implemented in 
some computational tools used for the analysis of the single-W^ process, nor 
the |g 2 * | scale are able to reproduce the effects due to the scales of eq. (14) and 
eq. (15). These two scales are in good agreement and both predict a lowering 
of the Born cross section of about 8-9%, almost independent of the cm. lep2 
energy and quark angular acceptance. This fact can be understood by looking 
at Fig. 5, where it is shown the single-W differential cross section with respect 
to the scales A± of eq. (14). On the left, A + exhibits a broad peak not far from 
M w , while, on the right, the other scale A_ peaks, as expected, at very small 
momentum transfer. 

Figure 6 shows the effects of the reweighting procedure of eq. (17) for the 
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Fig. 3. The effects of LL qed corrections to the cross section of the single- W process 
e + e~ — > e~vud for different choices of the energy scale in the electron/positron SF. 
The quark angular acceptance 0° < # j < 180° is considered. Left: absolute cross 
values as functions of the lep2 cm. energy. Right: relative difference between the 
qed corrected cross sections and the Born one, still as functions of the cm. energy. 
The marker • represents the Born cross section, o represents the correction given 
by Aj_ = s for both SF, o the correction given by the scales = \q^*\ for both SF, 
□ the correction given by the naive scales of eq. (15), A the correction given by the 
scales of eq. (14). The entries correspond to V* = 183, 189, 200 GeV. The markers 
are misplaced for better reading. 
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Fig. 4. The same as Fig. 3 for the quark angular acceptance | cos 9 U J\ < 0.95. 

evaluation of the QED running coupling constant. The marker A represents 
the relative difference between the integrated cross section computed in terms 
of ac F and the cross section computed in terms of cv(0), while the marker o is 
the relative difference between the integrated cross section computed in terms 
of «g f an d the cross section computed in terms of a(t). As can be seen, the 
rescaling from ac F to a(t) introduces a negative correction of about 5-6% in 
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Fig. 5. The differential cross sections of the single- W process e + e~ 
respect to the two scales of eq. (14) at yfs = 189 GeV. 
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Fig. 6. The effects of the rescaling of olqed from ac F to a(q^* = 0) (a) and 
a(q^*) (o) on the integrated cross section of the single- W process e + e~ — ► e~vud. 
o"o is the cross section computed in terms of clg f - The entries correspond to 
s/s = 183, 189, 200 GeV. 

the lep2 energy range. The difference between a and o, which is about 2-3%, 
is a measure of the running of olqed from <j^„ = to = t. 

As an illustrative example of the effect of anomalous couplings on single-VF 
differential distributions, in Fig. 7 the distribution of the qq invariant mass, 
around the peak of the VF-boson resonance, and the distribution of the angle 
of the quarks with the line of flight of the reconstructed VF-boson in the W- 
boson rest-frame are shown. The dashed lines correspond to the simulation 
as obtained by means of SWAP for the anomalous coupling Ak 7 = 0.1, while 
the solid lines represent the Standard Model prediction. The effect of the 
anomalous coupling Ak 7 at LEP2 energies is just an overall rescaling of the 
total cross section. Therefore the lep2 sensitivity to Ak 7 in single-VF events 
depends crucially on the accuracy of the theoretical evaluation of the total 
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Fig. 7. The single-TT differential cross sections with respect to the quark- antiquark 
invariant mass (left side), and with respect to the angle between a quark and the 
line of flight of the reconstructed VT-boson in the frame of rest of the TT-boson (right 
side). The dashed line represents the distribution in the presence of an anomalous 
gauge coupling Ak 7 = 0.1, while the solid line is the Standard Model prediction. 
The cm. energy is y/s = 189 GeV. 

cross section. 



7 Conclusions 

The process of single- W production in high-energy e + e~ collisions is relevant at 
lep2 for the determination of the non-abelian self-couplings of the W boson, 
and of primary importance at future Linear Colliders at the TeV scale, its 
cross section being dominant at very high energies with respect to other four- 
fermion processes. 

In order to give a contribution to the reduction of the theoretical uncertainty 
presently associated to the calculation of the single-!^ cross section, the issue 
of higher-order photonic corrections has been carefully investigated within the 
standard SF technique. Theoretical and phenomenological arguments for the 
choice of the energy scale entering the SF have been proposed. Two possible 
solutions for the scale of QED radiation have been obtained. The former has 
been derived by means of general arguments concerning the soft and collinear 
limit of the 0(a) corrections coming from the radiation of external legs. The 
latter, which can be considered as a naive ansatz, has been driven by thinking 
of the single-Vr process in terms of the Weizsacker- Williams approximation. 

Numerical calculations show that the typically adopted choice of the center- 
of-mass energy of the reaction, as radiation scale for the process, can lead to 
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over-estimate the radiative correction by a factor of 1.5, implying an under- 
estimate of the cross section of about 4%. Also the choice of fixing the scale 
to the momentum transfer t in the ecy* vertex for both the IS SF leads to an 
under-estimate of the photon correction of about 4%. The difference between 
the predictions given by the two set of scales of eq. (14) and eq. (15) is at 
the per mille level in the lep2 energy range. Therefore, the naive scales of 
eq. (15) provide a good ansatz for the energy scale of QED radiation in the 
single-Vr process, which could be simply implemented in MC tools for data 
analysis and further corroborated by the comparison with the results of other 
approaches. The method here described for the energy scale determination in 
the SF can be simply generalized to other four-fermion process dominated by 
non-annihilation channels, such as single- Z production. 

In order to provide adequate phenomenological predictions for precision ex- 
periments, also the running of the electromagnetic coupling constant has been 
accounted for in an effective way, i.e. by rescaling the differential cross section 
for the ratio of the electromagnetic coupling constant, evaluated at the typical 
scale of the process, to the same coupling evaluated from the input parame- 
ters according to tree-level relations. The effect of such rescaling amounts to 
a negative correction of about 5-6%, in agreement with recent findings [17], 
as far as the effect of ojqed is concerned. 

In the light of the experimental precision for the single-VT process, the correc- 
tions considered in the present paper are phenomenologically relevant. 

According to the theoretical approach described in the present paper, an orig- 
inal MC programme SWAP has been developed, including exact tree-level ma- 
trix elements with finite fermion masses effects, anomalous couplings, vacuum 
polarization and higher-order QED corrections. The code is available for ex- 
perimental analysis. 
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